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STRESS  CONCENTRATION  AROUND  BROKEN  FILAMENTS 
IN  A  FILAMENT-STIFFENED  SHEET  * 

By  W.  B.  Fichter 
Langley  Research  Center 

SUMMARY 

Stress  distributions  around  single  and  double  cuts  (groups  of  adjacent  broken  fila¬ 
ments)  in  an  idealized  composite  sheet  consisting  of  parallel  tension- carrying  filaments 
embedded  in  a  shear- carrying  matrix  are  analyzed  by  an  influence-function  technique. 
Static-load  concentration  factors  for  two  equal  collinear  cuts  are  obtained  and  are  com¬ 
pared  with  some  related  results  for  a  single  cut.  Dynamic-load  concentration  factors, 
corresponding  to  sudden  breaking  of  filaments,  also  are  obtained  for  two  of  the  simplest 
double-cut  cases.  In  addition,  matrix  shear  forces  are  investigated  for  some  single-cut 
and  double- cut  cases.  Finally,  loads  in  broken  filaments  are  calculated  for  some  single¬ 
cut  cases,  and  their  implications  in  some  current  failure  analyses  of  composite  materials 
are  briefly  discussed. 

The  interaction  between  two  collinear  cuts  is  essentially  local  and  is  roughly  con¬ 
fined  to  separation  distances  on  the  order  of  the  cut  length.  For  closely  spaced  cuts,  how¬ 
ever,  the  interaction  between  cuts  is  seen  to  be  pronounced.  For  a  single  cut,  an  expres¬ 
sion  analogous  to  the  filament  load  concentration  factor  is  obtained  for  the  maximuni  shear 
force  in  the  matrix  as  a  function  of  the  cut  length.  The  results  suggest  that  for  some 
combinations  of  constituent  mechanical  properties,  the  composite  might  be  more  suscep¬ 
tible  to  filament  tensile  failure  at  some  cut  lengths,  but  more  susceptible  to  matrix  shear 
failure  at  others. 


INTRODUCTION 

Much  effort  is  currently  being  directed  toward  the  development  of  lightweight  com¬ 
posite  materials  typically  composed  of  plastics  and  metals  in  which  are  embedded  high- 
strength  fibers.  An  extensive  survey  of  composite  materials  research  is  contained  in 


’This  paper  is  based  on  part  of  a  thesis  entitled  "Stress  Concentrations  in  Filament- 
Stiffened  Sheets"  submitted  in  partial  fulfillment  of  the  requirements  for  the  degree  of 
Doctor  of  Philosophy  in  Engineering  Mechanics,  North  Carolina  State  University,  August 
1969. 


reference  1.  In  studies  of  the  tensile  strength  of  fibrous  composites  (for  example, 
refSe  1  and  2),  it  has  been  observed  that  failure  can  be  initiated  by  local  cuts  or  breaks 
in  the  filaments.  A  theory  for  predicting  composite  tensile  strength  which  recognizes  the 
statistical  nature  of  filament  strength  as  well  as  local  cut  distribution  has  been  developed 
in  reference  2.  In  the  statistical  approach,  however,  the  stress  concentration  effects  of 
flaw  size  and  flaw  distribution  are  generally  ignored  or  treated  rather  approximately. 
Therefore,  an  analysis  in  which  stress  concentration  effects  are  taken  into  account  in 
more  detail  is  desirable,  and  the  present  paper  is  directed  toward  determining  some  of 
these  effects. 

Analytical  investigations  of  the  stress  fields  around  some  isolated  ''cutouts”  in  an 
infinite  flat  sheet  of  parallel  tension-carrying  filaments  embedded  in  a  shear-carrying 
matrix  are  contained  in  reference  3.  Included  in  these  investigations  is  an  analysis  which 
yields  the  filament  load  concentration  factor  for  a  single  straight  cut  across  an  arbitrary 
number  of  adjacent  filaments  (or,  equivalently,  a  group  of  adjacent  broken  filaments). 
More  recently,  this  analysis  has  been  extended  in  reference  4  to  problems  of  load  con¬ 
centration  in  composites  with  two-dimensional  arrays  of  parallel  filaments.  Another 
logical  extension  of  the  analysis  is  to  the  problem  of  two  cuts  in  the  same  filamentary 
sheet,  since  neighboring  flaws  in  a  filamentary  composite  can  be  expected  to  influence 
one  another. 

In  the  present  paper  an  analysis  is  given  of  the  stresses  in  a  filament  sheet  weak¬ 
ened  by  two  collinear  cuts  of  arbitrary  length.  For  the  particular  case  of  cuts  of  equal 
length,  filament  load  concentration  factors  are  computed  for  various  combinations  of  cut 
length  and  distance  between  cuts.  In  addition,  matrix  shear  forces  and  filament  tensile 
forces  are  computed  for  several  cases  involving  single  and  double  cuts.  For  the  case  of 
a  single  cut,  an  expression  analogous  to  the  filament  load  concentration  factor  is  obtained 
for  the  maximum  matrix  shear  force  as  a  function  of  cut  length.  Finally,  the  recovery  of 
load  by  broken  filaments  is  investigated  for  an  idealized  model  and  the  results,  as  they 
pertain  to  statistical  strength  predictions  such  as  those  of  reference  2,  are  briefly 
discussed. 


SYMBOLS 

d  filament  center-line  spacing 

df  filament  diameter 

EA  extensional  stiffness  of  a  filament 

Ej^  Weber  function 
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matrix  shear  modulus 


effective  thickness  of  matrix 


modified  Bessel  function  of  first  kind 


Bessel  function  of  first  kind 


i,k,m,n, 

q,r,s 


integers 


filament  load  concentration  factor 


modified  Struve  function 


mass  per  unit  length  associated  with  a  filament 

load  in  nth  filament  for  influence-function  solution  of  ref. 


dimensionless  load  in  nth  filament, 


force  applied  to  each  filament  at  infinity 


load  in  nth  filament 


maximum  dimensionless  matrix  shear  force 


dimensionless  matrix  shear  force  per  unit  length  between  nth  and  (n  +  l)th 
X., _ Sn  IeM 


filaments 


P  1/  Gh 


matrix  shear  force  per  unit  length  between  nth  and  (n  +  l)th  filaments 


dimensionless  displacement  of  nth  filament,  — 

P  V  d 


displacement  of  nth  filament 
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Vn 

X 

r 

T 

e 


displacement  of  nth  filament  for  iafluence-function  solution 
coordinate  parallel  to  filaments 
gamma  function 


dimensionless  coordinate  parallel  to  filaments, 


dim.ensionless  time, 


variable  of  integration 


GOVERNING  EQUATIONS 


Insofar  as  is  practical,  the  notation  employed  in  reference  3  is  retained  in  the 
present  analysis.  The  model  is  shown  in  figure  1(a),  along  v/ith  the  coordinate  system 
and  some  notation.  The  analytical  model  is  one  commonly  used  in  '^shear  lag’’  analyses, 
which  does  not  account  for  local  stress  variations  near  the  filament-matrix  interface.  It 
is  composed  of  parallel  tension- carrying  members  (filaments)  embedded  in  a  matrix  which 
carries  only  shear.  The  filaments  are  separated  by  a  constant  distance  d  and  are  num¬ 
bered  from  to  oo.  The  coordinate  along  the  filaments  is  x  and  the  displacement  of 
the  nth  filament  at  location  x  is  Un(x).  The  force  in  the  nth  filament,  taken  positive  in 
tension,  is  denoted  by  p^(x)  and  is  given  in  terms  of  u^^  by 

du„ 

p  (x)  =  EA  —  (1) 

dx  ^  ^ 

where  EA  is  the  extensional  stiffness  of  the  filament.  The  shear  force  per  unit  length 
between  the  nth  and  (n  +  l)th  filaments  is  defined  here  by  =  "^(^n+1  ”  Static 

equilibrium  of  an  element  of  the  nth  filament  then  requires 


In  fig-are  1(a),  i;ilaments  -2,  -1,  2,  and  3  are  shown  cut  at  x  =  0  for  illustrative 
purposes.  In  general,  for  two  collinear  cuts  through  q  and  s  filaments,  let 
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-(k  +  q)  ^  n  g  -(k  +  1)  and  r  +  lgn^r  +  s,  respectively,  denote  the  broken  filaments, 
the  two  cuts  being  separated  by  r  +  k  +  1  intact  filaments,  (See  fig.  1(b).)  Then  the 
boundary  conditions  at  x  =  0  are 

Uj^(O)  =  0  (-k  ^  n  ^  r;-(k  +  q  +  1)  ^  n;n  ^  r  +  s  +  l)  '1 

>  ( 

Pn(0)  =  0  (-(k  +  q)  g  n  g  -(k  +  l);r  +  lSnSr  +  s) 

For  X  large,  the  force  in  each  filament  approaches  the  uniform  applied  force,  denoted 
by  p;  that  is. 


Pj^(±°o)  =  p 


For  convenience,  let 


Pn  =  PPn 


^VEAGh 


/Gh  „ 

"n^ViAd^n 


Then  the  equilibrium  equations  and  boundary  conditions  become 


f +U^+l-  2Un  +  Un.i=0 


Un(0)  =  0  ^n  2  -(k  +  q  +  l);-k  g  n  ^  r;n  S  r  +  s  +  l) 

Pj^(O)  =  0  ^-(k  +  q)  5  n  s  -(k  +  l);r  +  lSngr  +  s) 


Pn(±°°)  =  1 


The  dimensionless  forces  and  displacements  are  related  by 


Sn  =  -  U„ 
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SOLUTIONS 


The  mixed  boundary  value  problem  defined  by  equations  (6),  (7),  and  (8)  can  be  solved 
by  use  of  the  influence-function  technique  employed  in  reference  3  to  solve  the  single-cut 

problem.  The  influence  functions  Vj^(^)  and  = ‘-^(^)  are,  respectively,  the  non- 

dimensional  displacement  and  force  in  the  nth  filament  when  the  filament  sheet  is  com¬ 
pletely  cut  along  ^  =  0  and  the  zeroth  filament  is  displaced  a  unit  amount  at  ^  =  0  while 
ail  other  filaments  are  held  fixed  at  ^  =  0.  In  terms  of  Vn(C)  ^n(^)  (^n(^) 

ref.  3),  the  dimensionless  force  and  displacement  in  the  nth  filament  are  given  by 


00 


In  reference  3,  Vj^(^)  has  been  found  to  be 

.  rTT 


pir  -2|sin  | 

V„(^)  =  -  \  cos  n0  e  d 

“  TT  Jq 

Application  of  boundary  conditions  (7)  to  equations  (9)  yields 

-(k+1)  r+s 


Pn(?)  =  l+  ^  Uj(0)  +  ^  Nn_j(|)Ui(0) 

i=-(k+q) 


(k+1) 


Un(l)  =  ?+  ^  Ui(0)  +  2,  Vi(?)Ui(0) 

i=-(k+q)  i=r+l 


values  of 

i,  and 

-(k+1) 

T - ^ 

r+s 

V — ' 

I  1^= 

„-i(0)Ui(0)+  2,  Nj,_i(0)  Ui(0) 

i=-(k+q) 

i=r+l 

(-(k  +  q)  g  n  s  -(k  +  l);r  +  1  g  n 

6 


which  expresses  the  condition  of  zero  load  on  the  ends  of  the  broken  filaments.  Equa¬ 
tions  (12)  constitute  a  set  of  q  +  s  linear  algebraic  equations  in  the  q  +  s  unknowns 
Uji(O)-  Their  solution  set  can  be  substituted  into  equations  (11)  to  obtain  expressions  for 
load  and  displacement  in  any  filament. 

However,  before  the  loads  and  displacements  can  be  calculated,  the  integral 
representation  of  the  influence  functions  Vjj(^)  and  Njj(^)  must  be  evaluated.  Only 

N  (0)  =  ^^(0)  has  been  evaluated  in  reference  3  because  the  computation  of  load  con- 
^  d^ 

centration  factors,  which  was  the  main  purpose  of  that  investigation,  does  not  require  the 
evaluation  of  the  influence  functions  for  nonzero  values  of  Equation  (10)  can  also  be 
written  as 


Vn(^)  =  ^ 


^  ^  .  -2^  Sin 

cos  2n0  e  ^ 


Integration  of  this  expression  gives 

Vn(^)  =  J2n(2i^)  '  iE2n(2i^) 

where  i  denotes  \PI,  J2jj  is  the  Bessel  function  of  the  first  kind  of  order  2n,  and 

Er,  is  the  Weber  function  of  order  2n.  (See  ref.  5.)  In  terms  of  functions  with  real 
2n 

arguments, 


Vn(^)  =  (-lfl2n(2?)-L2n(2^)-i>  (“D 


r(k.l>“-i 


r2nH-i-k 


(n  go) 


in  which  l2jj  is  the  modified  Bessel  function  of  the  first  kind,  L2n  is  the  modified 
Struve  function  (see  ref.  5),  and  the  finite  sum  must  be  taken  to  be  zero  for  n  =  0.  A 
similar  expression  for  V]fj(^)  can  be  derived  for  n  <  0;  however,  it  is  obvious  from 
equation  (13)  that  V.^d)  =  V^d).  Differentiation  of  equation  (14)  yields 


Nn(^)  =  ■  ^2n+l^^^^  L2n-l(2?)  +  1/2^(211  +  ^ 


n-1  (2n  -  2k  -  l)r  k  +  ik 

ly  (.i)k - V - ^ 

’’teO  r'(2„  +  l-k) 


l\.2n-2k-2 


(n  go) 


which  the  finite  sum  again  is  to  be  omitted  for  n  =  0.  From  equation  (15)  it  is  found 
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(16) 


Nn(0)  = 

77(411^  -  1) 

which  is  in  agreement  with  reference  3. 

With  substitution  of  the  appropriate  values  from  equation  (16),  equations  (12),  v/hich 
are  merely  linear  algebraic  equations,  can  be  solved  for  the  displacements  of  the  ends  of 
the  broken  filaments  11^(0).  Then  the  solutions  to  equations  (12),  along  with  the  influence 
functions  given  by  equations  (14)  and  (15),  can  be  substituted  into  equations  (11)  to  obtain 
the  load  and  displacement  at  any  point  in  the  filamentary  sheet.  However,  because  each 
cut  may  traverse  any  number  of  filaments  and  the  two  cuts  may  be  separated  by  any 
number  of  intact  filaments,  the  general  case  of  two  collinear  cuts  is  without  symmetry. 
Hence,  calculations  covering  reasonably  wide  variations  of  the  pertinent  parameters 
would  entail  considerable  computational  effort.  If  the  problem  is  simplified,  however,  by 
requiring  the  cuts  to  be  of  equal  length,  the  most  essential  features  of  the  two-cut  problem 
are  retained  and  the  computational  effort  is  greatly  reduced  through  consideration  of  the 
resulting  symmetries. 

In  what  follows,  then,  it  is  assumed  that  the  two  collinear  cuts  are  of  equal  length. 

In  the  analysis  of  this  reduced  problem,  one  of  two  cases  arises,  depending  upon  whether 
the  number  of  filaments  between  the  cuts  (henceforth  called  "interior  filaments")  is  even 
or  odd.  The  analysis  of  these  two  cases  is  presented  in  the  appendix,  pertinent  results 
being  drawn  from  there  as  appropriate. 

Filament  Loads 

In  the  case  of  a  single  cut  across  n  consecutive  filaments  (beginning  with  the 
zeroth  filament),  the  load  and  displacement  in  the  rth  filament  were  found  in  reference  3 
to  be 

n-1 

Pr(?)  =  1  +  ^  Nr_i(|)  Ui(0)  (17) 

i=0 


Ur(|)  =  I  +  ^  V^-i(?)  Ui(0)  (18) 

i=0 

where  the  Uj^(O)  terms  are  the  solutions  to  the  n  simultaneous  equations 

n-1 

0  =  1  +  ^  N^.i(O)  Ui(0)  (0  ^  r  s  n  -  1)  (19) 

i=0 
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For  the  case  of  two  equal-length  collinear  cuts  separated  by  any  number  of  intact  fila¬ 
ments,  expressions  for  the  load  and  displacement  in  any  filament  are  derived  in  the 
appendix. 

In  reference  3  the  load  concentration  factor  for  a  single  cut  across  n  filaments 
was  found  to  be 

n-1 

Kn  =  Pn(0)  =  1  ^  ^n-i^O)  Ui(0)  (20) 

i=0 

which  is  the  maximum  load  in  the  filament  just  ahead  of  the  cut.  For  double  cuts,  the  most 
highly  stressed  filaments  are  the  outermost  of  the  filaments  between  the  cuts  and,  again, 
their  maximum  load  occurs  at  ^  =  0.  In  the  appendix,  expressions  for  the  load  concen¬ 
tration  factors  have  been  found  for  cases  involving  odd  and  even  numbers  of  interior  fila¬ 
ments.  For  an  odd  number  2r  +  1  of  intact  filaments  between  two  cuts,  each  across  n 
filaments,  the  load  concentration  factor  is 

r+n 

Kn,2r+l  =  l+  J  |Nr+i(0)  +  N^.iCO^UiCO)  (21) 

i=r+l 

where  the  Ui(0)  terms  are  the  solutions  to  equations  (A4)  for  a  specific  value  of  n;  and 
for  an  even  number  2r  of  interior  filaments, 

r+n-1 

Kn,2r  =  l+  J  [Nr+i(0)  +  Nr-i-l(o3Ui(0)  (22) 

i=r 

where  the  Ui(0)  terms  are  calculated  from  equations  (AlO). 

Shear  Loads 

The  dimensionless  shear  force  per  unit  length  between  the  nth  and  (n  +  1)’-"  filaments 
is  given  by 

Sn(?)  =  Un+l(^)  -  Un(^)  (23) 

In  either  a  single-cut  or  double-cut  configuration,  the  greatest  matrix  shear  force  is 
between  a  broken  filament  and  an  unbroken  one.  In  the  case  of  a  double  cut,  the  unbroken 
filament  is  an  interior  one. 

To  find  the  matrix  shear  forces  in  a  single-cut  problem,  the  solutions  to  the  boundary 
load  conditions  (19)  are  substituted  into  equation  (18)  to  obtain  displacements,  and  those 
results  are  used  in  equation  (23).  In  a  double-cut  problem,  either  equations  (A4)  or  (AlO) 
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are  solved,  depending  on  the  number  of  interior  filaments,  and  the  solutions  are  substi¬ 
tuted  into  the  second  of  either  equation  (A3)  or  (A9),  respectively,  and  then  equation  (23) 
is  used. 


Calculations  of  matrix  shear  forces  have  been  made  for  some  single- cut  and  double¬ 
cut  cases.  The  calculations  are  restricted  in  each  case  to  the  matrix  area  between  two 
filaments  which  experiences  the  greatest  shear  stress.  The  solutions  are: 

For  single  cut  across  one  filament: 

So(^)  -  ||^  -  2[lo(2^)  -  Lo(20j  -  |[l_i(2^)  -  Ii(2?|  (24) 

where 

So(0)  =  -  I 

For  single  cut  across  two  filaments: 


-  ‘■o(2a]  +  1(2  +  -  Il(2?)J  - 


where 


Si(0)  = 


Stt 


For  single  cut  across  three  filaments: 


(25) 


S2(l)  =  -  f|/  4  +  m\[iQ(2|)  -  Lo(2C)]  +  1(20  + 
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+  ^YL_i(2|)-Ii(2?)1-i--I4  _6^ 


where 


S2(0)  = 


IStt 

32 


(26) 


For  two  cuts,  each  across  one  filament,  with  one  interior  filament: 

So<«  "  1^12(2  .  ^j[lo(2|)  -  Lo(25)].  |(l  .  ^)[l..i(2?)  -  Ii(2£j]  -  ±-  M 


where 


^  1577 


So(0)  =  5, 


(27) 
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For  two  cuts,  each  across  two  filaments,  with  one  interior  filament: 


So(^)  = 


SSGStt 

2(6512) 


where 


SqCo)  = 


55657r 

2(6512) 


and  where  the  bars  are  used  to  identify  double-cut  shear  forces.  A  comparison  of  equa¬ 
tions  (27)  and  (28)  will  indicate  how  rapidly  the  complexity  of  double-cut  shear  force  cal¬ 
culations  grows  with  cut  length. 


RESULTS  AND  DISCUSSION 


The  solutions  presented  in  the  preceding  section  have  been  employed  to  investigate 
various  phenomena  associated  with  the  breaking  of  filaments,  such  as  static  and  dynamic 
load  concentration  in  filaments,  matrix  shear  force  growth  with  cut  length,  as  well  as 
decay  of  matrix  shear  force  and  recovery  of  load  by  broken  filaments  as  functions  of  axial 
distance  from  the  cut. 


Load  Concentration  Factors 

Tensile  loads.-  In  reference  3  the  static  load  concentration  factor  for  a  single  cut 
across  n  consecutive  filaments  was  found  to  be 


3 


6  •  8  • 
5  •  7  • 


(2n  +  2) 
(2n  +  1) 


(n  =  1,  2,  3 


In  the  present  paper,  filament  load  concentration  factors  for  two  collinear  cuts  Kn^m 
each  across  n  filaments,  separated  by  m  intact  filaments,  have  been  calculated  for 
n  from  1  to  8  and  m  from  1  to  16.  The  results  are  presented  in  table  I  and  are  plotted 
in  figure  2,  along  with  the  single-cut  factor  Kn  for  purposes  of  comparison.  It  might 
be  noted  here  that  Kn  is  the  asymptotic  value  of  Kn^m  that  is,  for  large 

distances  between  the  cuts.  It  can  be  seen  in  figure  2  that  the  interaction  between  collin¬ 
ear  cuts  is  essentially  local,  being  confined  to  separation  distances  which  are  on  the  order 


of  the  cut  lengths  For  cuts  which  are  farther  apart,  the  results  do  not  differ  conspicuously 
from  their  asymptotic,  single-cut  values.  Thus,  when  two  cuts  are  separated  by  distances 
which  are  greater  than  their  lengths,  each  cut  can  be  considered  independently  with  little 
loss  of  accuracy. 

The  values  of  closely  spaced  cuts,  however,  increase  rapidly  when  m 

becomes  small.  Uence,  Kn^ni  small  is  of  particular  interest,  since  it  is  asso- 

cia.ted  with  states  of  severe  load  concentration.  In  order  to  examine  the  relafive  severity 
of  double  cuts  and  comparable  single  cuts,  double-cut  factors  for  closely  spaced  cuts  are 
plotted  in  figure  3  along  with  factors  for  single  cuts  of  comparable  total  length.  In  terms 
of  the  notation  employed  here,  i  and  o  compared  with  single-cut  factors 
denoted  by  K2p,  K2n+i,  and  K2ji+2-  These  particular  values  of  the  single-cut  factor 
are  employed  because  they  represent  either  the  same  total  number  2n  of  broken  fila¬ 
ments  or  the  number  2n  +1  or  2n  -f  2  which  would  be  broken  if  the  double  cuts  sepa¬ 
rated  by  one  or  two  intact  filaments,  characterized  by  Kj^  1  or  KjQ  2j  were  to  coalesce. 

It  can  be  seen  in  figure  3  that  the  load  concentration  factors  for  the  double  cut  with 
one  interior  filament  is  greater  than  those  for  a  single  cut  across  an  equal  number 

of  filaments  when  more  than  two  filaments  are  broken.  Also,  the  factors  for  the  double 
cut  with  two  interior  filaments  Kj-j^2  greater  than  those  for  a  single  cut  across  an 
equal  number  of  filaments  when  15  or  more  filaments  are  broken.  Thus,  it  appears  that 
the  longest  cut  does  not  necessarily  induce  the  most  severe  stress  state. 

This  result  suggests  that  a  composite  structural  design  criterion  based  on  residual 
strength  of  the  composite  in  the  presence  of  a  single  cut  of  prescribed  length  could  be 
unconservative,  even  though  only  cuts  of  lesser  length  were  present. 

In  reference  3  dynamic  load  concentration  factors  were  obtained  for  cases  in  which 
a  single  cut  is  suddenly  introduced  in  a  stretched  filament  sheet  a.nd,  in  a  separate  analy¬ 
sis,  an  apparent  upper  limit  of  1.27  v/as  found  for  the  dynamic  overshoot  (ratio  of  maxi¬ 
mum  dynamic  and  static  load  concentration  factors).  To  investigate  the  possibility  of  a 
departure  from  the  trend  for  a  single  cut,  dynamic  load  concentration  factors  have  been 
computed  for  two  of  the  simplest  double  cut  cases.  The  dynamic  analysis  was  carried 
out  by  inserting  the  appropriate  inertia  term  in  equation  (2)  and  utilizing  the  Laplace 
transform  theory  in  the  manner  of  reference  3.  The  results  are  presented  in  figure  4, 
v/here  som.e  single  cut  results  are  reproduced  for  purposes  of  comparison.  For  the  two 
double-cut  cases,  in  which  totals  of  two  and  four  filaments  were  suddenly  broken,  one 
remaining  intact  between  them,  the  dynamic  overshoots  were  found  to  be  1,22  and  1,23, 
respectively,  and  are  in  the  range  of  previously  obtained  single-cut  results.  The  fact 
that  additional  double-cut  calculations  ha.ve  not  been  made  because  of  their  rapidly 
increasing  complexity  precludes  the  drawing  of  sweeping  conclusions  concerning  the 
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overall  behavior  of  the  double-cut  dynamic  overshoot;  however,  the  present  limited  cal¬ 
culations  suggest  that  departure  from  the  range  of  single-cut  values  is  unlikely. 

Shear  loads.-  Although  the  present  composite  model  accounts  for  matrix  shear 
stress  in  only  an  approximate  manner,  it  is  believed  that  matrix  shear-force  calculations 
are  of  at  least  qualitative  interest.  A  result  which  was  not  obtained  in  reference  3,  but 
which  can  be  extracted  from  the  analysis  presented  there,  is  the  magnitude  of  the  most 
severe  shear  force  at  a  single  cut  as  a  function  of  the  number  of  cut  filaments,  which 
might  be  thought  of  as  a  shear-force  concentration  factor  Smax-  For  a  cut  which  starts 
at  the  zeroth  filament  and  severs  n  filaments,  the  most  severe  shear  force  is  given  by 

SmK  =  |S„(»)|  '  UoO) 

where  Sj^ax  defined  as  the  absolute  value  of  the  peak  shear  force,  and  where  Uo(0) 
is  the  nondimensional  displacement  of  the  end  of  the  zeroth  filament,  which  varies  with 
the  number  of  cut  filaments.  Solution  of  equations  (19),  (18),  and  (23)  for  the  first  six 
values  of  n  yields  the  following  values  of  Smax  ^^®o  eqs.  (24),  (25),  and  (26)): 


n 

C 

"^max 

1 

0.785 

2 

1.178 

3 

1.473 

4 

1.718 

5 

1.933 

6 

2.126 

Inspection  of  these  values  in  fractional  form  indicated  that  they  conformed  to  the 
expression 


4  o  (2n  -  D: 


(31) 


By  using  equations  (27)  and  (28),  a  few  results  were  also  obtained  for  the  shear  load 
concentration  factor  simplest  double-cut  configurations  in  which  a  sin¬ 

gle  intact  filament  separates  the  cuts.  Solutions  for  2  and  4  broken  filaments  yield 


— 

n 

®max 

2 

0.841 

4 

1.342 
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Unlike  tensile  load  concentration  factors  for  closely  spaced  cuts,  these  limited  results 
suggest  ths.t  the  shear  load  concentration  factors  for  double  cuts  are  less  than  those  for 
single  cuts  v;ith  the  same  number  of  cut  filaments. 


A  further  consideration  in  comparing  shear  and  tensile  load  concentration  factors 
is  their  relaiive  rates  of  increase  with  increasing  number  of  cut  filaments.  In  figure  5, 
the  shear  concentration  factor  for  single  cuts  is  compared  with  the  tensile  load 

concentration  factor  for  values  of  n  up  to  13.  The  shear  force  can  be  seen  to 
increase  more  rapidly  than  the  tensile  force,  particularly  for  small  values  of  n.  Both 
Kn  and  are  unbounded  as  n  approa.ches  infinity.  Hov/ever,  their  relaiive  mag¬ 

nitudes  for  large  n  can  be  determined  by  v;riting  equation  (29)  in  the  form 


T-  ■  2%’.(n  +  1): 

(2n  +  1): 

and  forming  the  ratio 

|Smarc  n(2n  +  l)[(2n)[|^ 

"2^"-l(n+  l)(ni^ 


(32) 


By  the  use  of  asymptotic  formulas  for  the  factorial  function,  it  is  found  that 

S 

max  <  ,  , 

,,  ^  1  (n  oo)  (uS) 

Equations  (32)  and  (33)  indicate  that  the  maximum  dimensionless  matrix  shear  force  and 
filament  tensile  force  are  of  comparable  magnitude  for  all  values  of  n.  They  are  approx¬ 
imately  equal  for  large  values  of  n,  although,  as  noted,  the  maximum  shear  force  initially 
increases  more  rapidly  than  the  tensile  force.  This  difference  in  slopes  for  small  n 
suggests  that  some  composite  materials  could  be  more  susceptible  to  tensile  failure  of 
the  filaments  when  only  small  flaws  are  present,  but  might  be  more  susceptible  to  matrix 
shear  failure  when  weakened  by  larger  flaws. 

The  longitudinal  variation  of  the  shear  loads  was  also  investigated.  In  figure  6,  the 
most  severe  matrix  shear  forces  for  single  cuts  across  one,  two,  and  three  filaments  (see 
eqs.  (24),  (25),  and  (26))  are  plotted  as  functions  of  the  axial  distance  from  the  cut.  Also 
shown  are  results  for  the  two  simplest  double-cut  configurations  (see  eqs,  (27)  and  (28)) 
in  v/hich  a  single  intact  filament  separates  the  cuts.  The  curves  are  normalized  with 
respect  to  their  peak  values,  which  can  be  obtained  from  the  table  of  amplitudes  in  fig¬ 
ure  6  along  with  the  pertinent  equation  from  equations  (5).  It  can  be  seen  that  the  peaic 
values  increase  with  cut  length.  The  normalized  curves  differ  from  one  another  in  their 
rates  of  decay  with  axial  distance  from  the  cut,  the  decay  generally  being  slower  for 
longer  cuts.  It  might  also  be  noted  that  the  double -cut  shear  loads  decay  somewhat 
more  rapidly  than  those  for  the  single  cuts. 
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Loads  in  Broken  Filaments 


For  a  single  cut  across  n  consecutive  filaments  (starting  with  the  zeroth  fila¬ 
ment),  the  load  in  the  ith  filament  is  given  by 

n-1 

Pi(^)  =  1  +  ^  Ni_k(l)  Uk(0) 

k=0 

where  the  n  broken  filaments  are  identified  by  i  =  0,  1,. .  .  n  -  1,  and  each  11^(0)  is 
one  of  a  unique  set  associated  with  a  specific  value  of  n.  The  11^(0)  values  are  obtained 
by  solution  of  equations  (19). 

Loads  in  the  broken  filaments  have  been  calculated  for  a  single  cut  across  one,  two, 
and  three  filaments.  The  results  are  plotted  in  figure  7,  where  it  can  be  seen  that  the 
length  over  which  a  broken  filament  carries  less  than  a  particular  fraction  of  its  load  can 
vary  considerably  with  the  number  of  broken  filaments.  The  90-percent  load-recovery 
line  has  been  added  for  later  reference.  With  slightly  more  labor,  similar  results  can 
be  obtained  for  various  double-cut  configurations;  however,  results  for  double  cuts  are 
not  expected  to  exhibit  significantly  different  trends.  For  this  reason,  and  because  illus¬ 
tration  of  the  effect  of  cut  length  on  filament  load  recovery  is  the  main  purpose  of  the  cal- 
cvilations,  loads  in  broken  filaments  have  not  been  calculated  for  double  cuts. 

Ineffective  Length  Calculation 

In  reference  2,  the  ineffective  length  of  a  filament  is  defined  as  that  portion  of  a 
broken  filament  which  supports  less  than  90  percent  of  the  load  applied  at  infinity.  Studies 
cited  in  references  1  and  2  suggest  that  the  ineffective  length,  based  on  a  load-recovery 
fraction  of  90  percent,  of  a  single  broken  filament  embedded  in  a  matrix  can  range  from 
less  than  one  up  to  several  hundred  filament  diameters,  depending  on  the  geometry  and  the 
mechanical  properties  of  the  constituents.  Calculations  have  been  made  in  reference  2, 
for  example,  of  ineffective  length  as  a  fmiction  of  filament  volume  fraction  vj,  and  E/G, 
the  ratio  of  filament  Young's  modulus  to  matrix  shear  modulus.  Results  are  presented  in 
the  form  of  a  family  of  curves,  each  member  of  which  corresponds  to  a  particular  filament 
volume  fraction. 

In  the  case  of  a  composite  containing  a  cut,  an  additional  parameter  which  would  be 
expected  to  influence  the  ineffective  length  is  the  number  of  broken  filaments.  With  the 
aid  of  an  idealized  model,  results  from  the  present  study  can  be  used  to  obtain  an  indica¬ 
tion  of  the  influence  of  this  parameter.  A  typical  cross  section  of  this  model  is  shown  in 
the  sketch.  For  the  present  ineffective  length  calculations,  it  is  assumed  that  the  thick¬ 
ness  of  the  sheet  is  equal  to  the  filament  diameter,  and  that  in  the  matrix  between  any  two 
adjacent  filaments  the  shear  stress  is  constant  at  a  given  axial  location.  Then  for  this 
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model  the  filament  volume  fraction  is  given  by  Vf  =  -7-^  where  df  is  the  diameter  of  a 

^  4d  ^ 

filament  and  d  is  the  filament  spacing. 

In  terms  of  the  pertinent  mechanical  and  geometrical  parameters,  the  axial  distance 
along  a  filament  is  given  by 


which  for  the  idealized  model  becomes 


In  order  to  obtain  estimates  of  ineffective  length  changes  due  to  the  breaking  of  additional 
filaments,  it  is  necessary  to  determine  from  figure  7  the  values  of  denoted  by  at 
which  the  broken  filaments  of  interest  have  recovered  a  stipulated  fraction  of  their  far- 
field  load.  These  values  of  |,  along  with  appropriate  values  of  E/G  and  Vj,  are  then 
substituted  into  equation  (34)  to  obtain  values  of  the  nondirnensional  ineffective  length. 


denoted  by  xy^d^ 


This  procedure  ha.s  been  ca.rried  out  for  two  values  of  V£  and  a  range  of  values  of 
E/G,  the  filament  being  considered  ineffective  wherever  it  supports  less  than  90  percent 
of  its  far-field  load.  The  cases  of  cuts  across  one,  two,  and  three  adjacent  filaments  have 
been  considered.  The  results  are  presented  in  figure  8.  In  the  case  of  the  cut  across 
three  filaments,  the  90-percent  recovery  figure  pertains  to  either  of  the  two  outer 
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filaments  in  the  cut  rather  than  to  the  middle  filament  which  would  have  yielded  even 
greater  values  of  ineffective  length. 

The  mathematical  model  employed  in  the  analysis  of  reference  2  consists  of  a  sin¬ 
gle  filament  encased  in  a  thin  layer  of  shear-carrying  material  (binder)  which,  in  turn, 
is  embedded  in  an  infinite  body  to  which  are  assigned  the  average  stiffness  properties  of 
the  composite  material.  Hence,  interaction  between  neighboring  filaments  is  ignored.  In 
addition,  shear  stresses  in  the  average  material  are  assumed  to  decay  in  a  negligible  dis¬ 
tance  from  its  interface  with  the  binder  material.  Because  of  these  basic  differences 
with  the  present  model,  the  results  of  reference  2  are  not  plotted  in  figure  8.  It  should 
be  noted,  however,  that  the  results  of  reference  2  yield  a  single  curve  for  each  filament 
volume  fraction,  whereas  the  present  results  yield  a  family  of  curves,  each  member  of 
which  corresponds  to  a  particular  number  of  broken  filaments.  In  addition,  the  analysis 
of  reference  2  predicts  a  more  rapid  recovery  of  load  by  broken  filaments  than  does  the 
present  analysis.  As  can  be  seen  in  figure  8,  large  changes  in  filament  ineffective  length 
can  result  from  varying  the  number  of  broken  filaments.  The  results  suggest  that  this 
effect  should  perhaps  be  accounted  for  in  statistical  studies  of  the  strength  of  filamentary 
composites. 


CONCLUDING  REMARKS 

An  analysis  has  been  made  of  the  loads  and  deformations  in  a  filament- stiffened 
sheet  which  has  been  weakened  by  two  collinear  cuts.  Also,  some  results  for  the  case  of 
a  single  cut,  in  addition  to  those  previously  reported,  have  been  obtained.  It  has  been 
found  that  significant  interaction  between  collinear  cuts  is  largely  restricted  to  cases  in 
which  the  distance  between  cuts  is  no  greater  than  the  cut  length.  It  is  seen  that  two 
closely  spaced  cuts  can  cause  greater  stress  concentration  than  a  single  cut  across  a 
comparable  total  number  of  filaments.  Limited  calculations  of  dynamic  stress  concen¬ 
tration  factors  for  suddenly  introduced  collinear  cuts  support  the  earlier  conclusion  that 
dynamic  effects  are  not  of  great  importance  in  filamentary  composites  of  the  type 
investigated. 

Calculations  of  loads  in  broken  filaments  show  significant  changes  in  ineffective  fila¬ 
ment  length  with  cut  length;  thus,  statistical  strength  analyses  of  composite  materials 
should  consider  the  incorporation  of  a  flaw-size  parameter. 

Calculations  for  single  cuts  of  various  lengths  show  that  maximum  shear  forces  grow 
more  rapidly  with  cut  length  than  maximum  tensile  forces.  Such  calculations  might  be 
useful  in  determining  whether  a  composite  material  is  more  susceptible  to  shear  failure 
or  tensile  failure. 
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The  present  analysis  is  based  on  linear  small-deflection  theory  of  a  single  sheet, 
whereas  filamentary  composites  usually  are  many  filaments  thick  and  are  subject  to  var¬ 
ious  nonlinear  effects  such  as  plastic  deformations,  large  deflections,  and  straightening  of 
the  filaments.  Although  these  effects  would  be  expected  to  exert  some  perturbing  influ¬ 
ence,  the  present  results  would  be  expected  to  remain  qualitatively  descriptive  of  load 
concentration  effects  in  filamentary  composites  of  the  type  investigated, 

Langley  Research  Center, 

National  Aeronautics  and  Space  Administration, 

Langley  Station,  Hampton,  Va.,  June  30,  1969. 
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APPENDIX 


STRESS  CONCENTRATION  FACTORS  FOR  TWO 
EQUAL- LENGTH  COLLINEAR  CUTS 


In  this  appendix  attention  is  restricted  to  problems  in  which  the  two  collinear  cuts 
are  of  equal  length.  This  restriction  is  effected  by  setting  q  equal  to  s  in  equa¬ 
tions  (11)  and  (12);  this  procedure  yields 


and 


-(k+1) 

r+s 

V - V 

i  '"n- 

id)  Ui(0)  + 

I 

i=(k+s) 

i=r+l 

-(k+1) 

r+s 

I 

.i(^)  Ui(0) 

i=-(k+s) 

i=r-hl 

-(k+1) 


r+s 


0  =  1+  ^  Nn_i(0)  Ui(0)  +  ^  Nn_i(0)Ui(0) 

i=-(k+s) 


i=r+l 


(Al) 


(-(k  +  s)  =  n  s  -(k  +  l);r  +  lgn^r  +  s)  (A2) 

One  of  two  cases  arises,  depending  on  whether  the  number  of  intact  filaments  between  the 
cuts  (interior  filaments)  is  odd  or  even.  The  case  of  an  odd  number  of  interior  filaments 
is  treated  first. 


Odd  Number  of  Interior  Filaments 

When  the  two  cuts  are  of  equal  length,  the  displacements  are  symmetric  with  respect 
to  the  line  which  equally  divides  the  group  of  interior  filaments.  In  this  case  it  is  con¬ 
venient  to  take  the  zeroth  filament  as  the  line  of  symmetry  so  that  k  =  r  in  equations  (Al) 
and  (A2)  and  the  total  number  of  interior  filaments  is  2r  +  1.  With  Uj(0)  =  U_j(0) 
because  of  symmetry,  equations  (Al)  and  (A2)  become  after  some  manipulation 
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and 


r+s  ■>! 

Pna)  =  l+  ^  [Nn+i(^)  +  N,^-i(^)]Ui(0) 

i=r-{-l 

r+s 

Un(l)  =  l+  ^  |yn+i(^)  +  Vn-i{?3Ui(0) 

i=r+l 


(A3) 


r+s 

0  =  1+  [Nn+i(0)  +  Nn_i(0)]Ui(0)  (r  +  1  g  n  ^  r  +  s)  (A4) 

i=r+l 

The  maximum  load  will  occur  adjacent  to  the  cuts  =  0)  in  the  two  outermost  interior 
filaments  (n  =  ±r).  Hence,  from  the  first  of  equations  (A3),  the  stress  concentration  fac¬ 
tor  for  two  cuts  across  s  filaments  separated  by  2r  +  1  filaments  is 

r+s 

Ks,2r+1  =  Pr(0)  =  1  +  ^  j^r+i^®)  +  Nr_i(0^Ui(0)  (A5) 

i=r+l 

where  the  Ui(0)  terms  are  calculated  from  equation  (A4)  for  specified  values  of  r 
and  s. 


Even  Number  of  Interior  Filaments 

In  this  case,  the  interior  filaments  are  assumed  to  be  identified  by  -k  g  n  2  r, 
which  with  r  =  k  -  1  becomes  -k  £  n  g  k  -  1.  Equations  (Al)  and  (A2)  take  the  form 


-(k+1) 

k+s-1 

Pn(l)  =  1  +  2 

Nn-i(l)Ui(0)+  ^ 

Nn-i(l)  Ui(0) 

i=-(k+s) 

i=k 

-(k+1) 

_ — , 

k+s-1 

+ 

•U/' 

1! 

Vn-i(?)  Ui(0)  + 

Vn-i(l)  Ui(0) 

i=-(k+s) 

i=k 

and 
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-(k+1)  k+s-1 

0  =  1+  ^  Nn-i(O)  Ui(0)  +  ^  Nn-i(O)  Ui{0) 

i=-(k+s) 

(-(k  +  s)  ^  n  ^  -(k  +  l);k  i  n  s  k  +  s  -  1))  (A7) 

For  this  configuration,  the  axial  line  of  symmetry  is  located  between  the  filaments  num¬ 
bered  0  and  -1,  which  means  that 

U_.(0)  =  Ui_i(0)  (A8) 

After  substitution  from  equation  (A8)  and  some  manipulation,  equations  (A6)  and  (A7) 
become 


r+s-1 

i=r 

r+s-1 

Una)  =  ?+  J  f'n+i+l(^)+Vn_ia)]Ui(0) 

i=r 


(A9) 


r+s-1 


0  =  1+  ^  jNjj^.j^j(0)  +  Nj^_|(0^U^(0) 


(r£nSr  +  s-l)  (AlO) 


In  this  case  the  maximum  load  occurs  for  n  =  -r,  r  -  1.  Then  with  n  -  r  -  1  in  the 
first  of  equations  (A9),  the  stress  concentration  factor  for  two  cuts  across  s  filaments, 
separated  by  2r  filaments,  is 


±  TO—  -L 

^(0)  =  1  +  ^  |Nr+i(0)  +  Nj._i.i(0gUi(0) 


(All) 


where  the  11^(0)  terms  are  calculated  from  equations  (AlO)  for  specified  values  of  r 
and  s. 
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Filament 
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(b)  Colli  near  cuts  across  q  and  s  filaments,  separated  by  r+k+1  intact  filaments. 
Figure  L-  Concluded. 


Figure  3.-  Comparison  of  filament  load  concentration  factors  for  single  and  double  cuts. 


Double  cut 


Figure  4.-  Dynamic  load  concentration  factors. 
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Figure  Ir  Loads  in  broken  filaments  near  a  single  cut. 


